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ABSTRACT 



We have made measurements of the Sunyaev-Zel'dovich (SZ) effect in six 
galaxy clusters at z > 0.2 using the Sunyaev-Zel'dovich Infrared Experiment 
(SuZIE II) in three frequency bands between 150 and 350 GHz. Simultaneous 
multi-frequency measurements have been used to distinguish between thermal 
and kinematic components of the SZ effect, and to significantly reduce the effects 
of variations in atmospheric emission which can otherwise dominate the noise. 
We have set limits to the peculiar velocities of each cluster with respect to the 
Hubble flow, and have used the cluster sample to set a 95% confidence limit 
of < 1410kms _1 to the bulk flow of the intermediate-redshift universe in the 
direction of the CMB dipole. This is the first time that SZ measurements have 
been used to constrain bulk flows. We show that systematic uncertainties in 
peculiar velocity determinations from the SZ effect are likely to be dominated by 
submillimeter point sources and we discuss the level of this contamination. 
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1. Introduction 

The spectral distortion to the Cosmic Microwave Background radiation (CMB) caused 
by the Compton scattering of CMB photons by the hot gas in the potential wells of galaxy 
clusters, known as the Sunyaev-Zel'dovich (SZ) effect, is now relatively straightforward to 
detect and has now been measured in more than 50 sources (see Carlstrom, Holder & Reese 
2002, for a review). Single-frequency observations of the SZ effect can be used to determine 
the Hubble constant (Myers et al. 1997; Holzapfel et al. 1997a; Pointecouteau et al. 1999; 
Reese et al. 2000; Pointecouteau et al. 2001; Jones et al. 2001; De Petris et al. 2002; Reese 
et al. 2002) and to measure the baryon fraction in clusters (Myers et al. 1997; Grego et al. 
2001). 

The spectrum of the SZ effect is also an important source of information. It can be 
approximated by the sum of two components (see Figure 1) with the strongest being the 
thermal SZ effect that is caused by the random thermal motions of the scattering electrons 
(Sunyaev & Zel'dovich 1972). The kinematic SZ effect, due to the peculiar velocity of 
the intracluster (IC) gas with respect to the Hubble flow (Sunyeav & Zel'dovich 1980), is 
expected to be much weaker if peculiar velocities are less than 1000 km s -1 , as favored by 
current models (Gramann et al. 1995; Sheth & Diaferio 2001; Suhhonenko & Gramann 2002). 
The thermal SZ effect has a distinct spectral signature, appearing as a decrement in intensity 
of the CMB below a frequency of ~ 217 GHz, and an increment at higher frequencies (the 
exact frequency at which the thermal effect is zero depends on the temperature of the IC 
gas, as discussed by Rephaeli 1995). The kinematic effect appears as a decrement at all 
frequencies for a cluster that is receding with respect to the Hubble flow, and an increment 
at all frequencies for a cluster that is approaching. Measurements that span the null of 
the thermal effect are able to separate the two effects, allowing the determination of the 
cluster peculiar velocity (Holzapfel et al. 1997b; Mauskopf et al. 2000; LaRoque et al. 2002). 
Additionally, SZ spectral measurements can, in principle, be used to determine the cluster 
gas temperature independently of X-ray measurements (Pointecouteau, Giard & Barret 1998; 
Hansen, Pastor & Semikov 2002), the CMB temperature as a function of redshift (Rephaeli 
1980; Battistelli et al. 2002) and also to search for populations of non-thermal electrons 
(Shimon & Rephaeli 2002). 

Peculiar velocities probe large-scale density fluctuations and allow the distribution of 
matter to be determined directly without assumptions about the relationship between light 
and mass. Measurements of a large sample of peculiar velocities can be used to probe Q m 
independently of the properties of dark energy (Peel & Knox 2002), and can, in principle, 
be used to reconstruct modes of the gravitational potential (Dore, Knox & Peel 2002). The 
local (z< 0.05) peculiar velocity field has already been measured and has been used to place 
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tight constraints on Vt m (Branchini et al. 2000; Courteau, Strauss & Willick 2000; Courteau 
& Dekel 2001; Bridle et al. 2001). However, the techniques used cannot be easily extended to 
higher redshifts because optical methods of distance determination have errors that increase 
linearly with distance. SZ spectral measurements allow peculiar velocities to be determined 
independently of the extragalactic distance ladder. 

In this paper we describe the measurements of the SZ spectrum of 6 galaxy clusters 
made with the Sunyaev-Zel'dovich Infrared Experiment (SuZIE). The SuZIE II receiver 
makes simultaneous measurements of the SZ effect in three frequency bands, centered at 
145 GHz, 221 GHz, and 355 GHz (or 270 GHz), spanning the null of the SZ thermal effect. 
We use these measurements to set limits on the peculiar velocity of each cluster, and to 
determine the cluster optical depth. The layout of this paper is as follows: in §2 we define 
our notation for the SZ effect; in §3 and §4 we describe the observations and data analysis. In 
§5 we consider statistical and systematic sources of uncertainty; in §6 we estimate the level 
of astrophysical confusion in our results and in §7 we determine the limits that the SuZIE 
peculiar velocity sample can be used to set on bulk flows, and discuss future prospects for 
this technique. 

2. The Sunyaev-Zel'dovich Effect 

2.1. The Thermal Effect 

We express the CMB intensity difference caused by a distribution of high energy elec- 
trons, n e , along the line of sight as (Rephaeli 1995): 

A/ T = J ^-[$( a ;,r e )-l]r (1) 
e x — 1 

where x = hu/kT , I = 2(kT ) 2 / (he) 2 , T is the temperature of the CMB, r — J n e a T dl 
is the optical depth of the cluster to Thomson scattering, and &(x,T e ) is an integral over 
electron velocities and scattering directions that is specified in Rephaeli (1995). In the limit 
of non-relativistic electrons, this reduces to the familiar non-relativistic form for the thermal 
SZ effect: 

Mx,T e )-l] = -^- 
e x — 1 

Following Holzapfel (1997b), we define: 

*(*,r e ) = ^[$(s,T e )-l] (3) 
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There also exist other analytic and numerical expressions for equation (1) (see Challinor 
& Lasenby 1998; Itoh et al. 1998; Dolgov et al. 2001) based on a relativistic extension of 
the Kompaneets equation (Kompaneets 1957). These expressions are in excellent agreement 
with equation (1). 

We define Comptonization as y = r x (kT e /m e c 2 ). It is a useful quantity because it 
represents a frequency independent measure of the magnitude of the SZ effect in a cluster 
that, unlike AI T , allows direct comparisons with other experiments. 



2.2. Kinematic SZ Effect 

The change in intensity of the CMB due to the non-relativistic kinematic SZ effect is: 

AI K = -J x x J n e a T ^ • dl (4) 

where v p is the bulk velocity of the IC gas relative to the CMB rest frame, and c is the 
speed of light in units of kms -1 . This functional form for the kinematic SZ effect has the 
same spectral shape as primary CMB anisotropy anisotropies, which represent a source of 
confusion to measurements of the kinematic effect. 

An analytic expression for the relativistic kinematic SZ effect has been calculated by 
Nozawa et al. (1998) as a power series expansion of 9 e = kT e /mc 2 and (3 = v p /c, where 
v p = v p ■ 1 is the radial component to the peculiar velocity. They found the relativistic 
corrections to the intensity to be on the order of +8% for a cluster with electron temperature 
kT e = 10 keV and v p = 1000 km s _1 . Although this is a relatively small correction to our 
final results, we use their calculation in this paper. We express the spectral shift due to the 
kinematic SZ effect as: 

AI K = -I xtx x h(x, T e ) (5) 

where h(x, T e ) is given by: 

h(x, T e ) = ^ x [1 + 9 e C 1 (x) + 9 2 e C 2 (x)} (6) 

and 9 e = kT e /m e c 2 . This expression includes terms up to 0(j39 2 e ). The quantities C±(x) 
and C 2 (x) are fully specified in Nozawa et al. (1998), who have also calculated corrections 
to equation (6) up to 0(j3 2 ). They find the correction from these higher order terms to be 
+0.2% for a cluster with kT e = 10 keV and v p = 1000 km s _1 , at a level far below the 
sensitivity of our observations. Therefore it can be safely ignored. 
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3. S-Z Observations 
3.1. Instrument 

In this paper we report measurements of the Sunyaev-Zeldovich effect made with the 
second generation Sunyaev-Zeldovich Infrared Experiment receiver (SuZIE II) at the Caltech 
Submillimeter Observatory (CSO) located on Mauna Kea. The SuZIE II receiver, described 
in Mauskopf et al. (2000) and associated references, measures the SZ effect simultaneously 
in three frequency bands. Two of the bands are centered at 145 and 221 GHz. The third was 
originally located at 273 GHz but has since been moved to 355 GHz to improve the degree 
to which correlated atmospheric noise can be removed from the data (Mauskopf et al. 2000). 
Between November 1996 and November 1997, the SuZIE II optics were changed slightly, 
altering both the beam size and chop throw. For a summary of the SuZIE II pass-bands 
and beam sizes in both configurations, see Table 1. The observations discussed in this paper 
include data in both configurations. 

The SuZIE II instrument consists of a 2 x 2 arrangement of 3-color photometers that 
observe the sky simultaneously in each frequency band. Each frequency is detected with 
a 300mK NTD Ge bolometer (Mauskopf et al. 1997). Light is coupled to the photometers 
through Winston horns which over-illuminate a 1.6K Lyot stop placed at the image of the 
primary mirror formed by a warm tertiary mirror. Each photometer defines a ~ 1'5 FWHM 
beam, with each row separated by ~ 2'.3 and each column by ~ 5' on the sky (see Figure 2). 
The beam size was chosen to correspond to typical cluster sizes at intermediate redshift 
(0.15 < z < 0.8). Bolometers in the same row that are sensitive to the same frequency are 
differenced electronically to give an effective chop throw on the sky of 5'. This reduces the 
level of common-mode atmospheric emission as well as common mode bolometer tempera- 
ture, and amplifier gain, fluctuations. This differencing strategy has been discussed in detail 
by Holzapfel et al. (1997a) and Mauskopf et al. (2000). 

3.2. Observation Strategy 

Observations of six clusters were made with SuZIE II over the course of several observing 
runs between April 1996 and November 2000 and are summarized in Table 2. As shown in 
Figure 2, SuZIE II operates in a drift scanning mode, where the telescope is pointed ahead 
of the source and then parked. The earth's rotation then causes the source to drift through 
the array pixels. Before each scan the dewar is rotated so that the rows of the array lie 
along lines of constant declination. Each scan last two minutes, or 30' in right ascension, 
during which time the telescope maintains a fixed altitude and azimuth. After a scan is 
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complete, the telescope reacquires the source and the scan is then repeated. Keeping the 
telescope fixed during an observation prevents slow drifts from changes in ground-spillover 
from contaminating the data. From scan to scan the initial offset of the telescope from 
the source (referred to in Figure 2 as RAOO and RAOl) is alternated between 12' and 18', 
allowing a systematic check for an instrumental baseline and a check for any time dependant 
signals. During the observations presented here, the array was positioned so that one row 
passed over the center of each cluster, as specified in Table 2. 



3.3. The Cluster Sample 

We selected bright, known X-ray clusters from the ROS 'AT X-Ray Brightest Abell Clus- 
ters (Ebeling et al. 1996a,b) and Brightest Cluster Samples (Ebeling et al. 1998) and the 
Einstein Observatory Extended Medium Sensitivity Survey (Gioia et al. 1990). The obser- 
vations of A1835 in April 1996 were previously analyzed by Mauskopf et al. (2000). We are 
using a variation of the analysis method used by Mauskopf and choose to re-analyze this 
data to maintain consistency between this cluster and the other data sets. We also use a 
slightly different X-ray model for this cluster, based on a joint analysis of radio SZ and X-ray 
data by Reese et al. (2002). Four of our clusters have been observed at 30 GHz by Reese 
et al. (2002). A comparison of our results with their measurements will be the subject of a 
separate paper. 



3.4. Calibration 

We use Mars, Uranus and Saturn for absolute calibration and to measure the beam 
shape of our instrument. The expected intensity of a planetary calibrator is: 

f2k(^T^ m (u)f k (u)du 

with T p i an (z/) being the Rayleigh- Jeans (RJ) temperature of the planet, and fk{v) the trans- 
mission function of channel k, whose measurement is described later in this section. We 
correct for transmission of the atmosphere by measuring the opacity using a 225 GHz tip- 
ping tau-meter located at the CSO. This value is converted to the opacity in each of our 
frequency bands by calculating a scaling factor a fe which is measured from sky dips during 
stable atmospheric conditions. For our frequency bands at 145, 221, 273, and 355 GHz we 
find a = 0.8, 1.0, 2.7, and 5.8. From drift scans of the planet we measure the voltage, Vpe^ 
that is proportional to the intensity of the source. We then find our responsivity to a celestial 
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source is: 



plan^ 'plan 



X e 



a<r/ cos 6» Cal > 



R = 



IV J 



(8) 



peak 



where f2 p i an is the angular size of the planet, and < rj costal > is averaged over the length 
of the observation, typically less than 20 minutes. We observe at least one calibration source 
every night. 

The data are then calibrated by multiplying the signals by a factor of R x e a<T ^ cosd sz> _ 
We correct for the transmission of the atmosphere during each cluster observation using 
the same method as for the calibrator observations. Each cluster scan is multiplied by 
e a<T l cos 6>sz> , where < rj cos#sz > is averaged over the length of that night's observation of 
the cluster, typically less than three hours. We average the atmospheric transmission over 
the observation period to reduce the noise associated with the CSO tau-meter measurement 
system (Archibald et al. 2002). To determine whether real changes in r over this time 
period could affect our results, we use the maximum variation in < rj cos#sz > over a single 
observation, and estimate that ignoring this change contributes a ±2% uncertainty in our 
overall calibration. 

The uncertainty of J p i an is dominated by uncertainty in the measurement of T p i an (z/). 
Measurements of RJ temperatures at millimeter wavelengths exist for Uranus (Griffin & 
Orton 1993), Saturn and Mars (Goldin et al. 1997). Griffin & Orton model their measured 
Uranian temperature spectrum, Tu raims (i>) , with a third order polynomial fit to the logarithm 
of wavelength. They report a 6% uncertainty in the brightness of Uranus. Goldin et al. 
report RJ temperatures of Mars and Saturn in four frequency bands centered between 172 
and 675 GHz. From these measurements we fit a second order polynomial in frequency to 
model Tg aturn (z/) and a second order polynomial in the logarithm of wavelength to model 
^Mars^)- Goldin et al. report a ±10K uncertainty to the RJ temperature of Mars due to 
uncertainty from their Martian atmospheric model, which translates to a 5% uncertainty in 
the brightness. They then use their Martian calibration to cross-calibrate their measurements 
of Saturn. Not including their Martian calibration error, they report a ~2% uncertainty 
to Saturn's RJ temperature. Adding the Martian calibration error in quadrature yields a 
total 5% uncertainty to the brightness of Saturn. The rings of Saturn have an effect on its 
millimeter wavelength emission which is hard to quantify. By cross-calibrating our Saturn 
measurements with SuZIE II observations of Mars and Uranus made on the same night, we 
have determined that ring angles between ±9° have a negligible effect on the total emission 
from Saturn. At higher ring angles we use Saturn only as a secondary calibrator. Since we 
use a combination of all three planets to calibrate our data, we estimate an overall ±6% 
uncertainty to T p i ane t. 

Bolometers have a responsivity that can change with the amount of power loading on 
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the detector; such non-linearities can potentially affect the results of calibration on a bright 
planet and the response of the detectors during the course of a night. We have used laboratory 
measurements to determine the dependence of responsivity on optical power loading. We 
estimate the variation in our loading from analysis of sky-dips taken at the telescope, and the 
calculated power received from Saturn, which is the brightest calibrator that we use. Over 
this range of loading conditions the maximum change in detector response is ~ 7.0%, 8.0%, 
and 3.5% in our 145, 221, and 355 GHz frequency bands respectively. Since the responsivity 
change will be smaller for the majority of our observations, we assign a 6% uncertainty in 
the calibration error budget for this effect. 

Mauskopf et al. (2000) has found that the SuZIE II beam shapes have a systematic 
dependence on the rotation angle of the dewar, which affects the overall calibration of the 
instrument. Based on these measurements we assign a ±5% calibration uncertainty from 
this effect. Further uncertainty to the calibration arises from our measurement of the spec- 
tral response, fk(v), which affects both the intensity that we measure from the planetary 
calibrators and the SZ intensity. The spectral response of each SuZIE II channel was mea- 
sured with a Michelson Fourier Transform Spectrometer (FTS). We then use the scatter of 
the measurements of the four bolometers that measure the same frequency to estimate the 
uncertainty in the spectral calibration at that frequency. We estimate this uncertainty to be 
<1% of the overall calibration. 

The calibration uncertainties are summarized in Table 3. Adding all of these sources in 
quadrature, we estimate the total calibration uncertainty of SuZIE II in each of its spectral 
bands to be ±10%. 



3.5. Definition of the Data Set and Raw Data Processing 

We now define some notation. Each photometer contains three bolometers each observ- 
ing at a different frequency. During a scan, the six difference signals that correspond to the 
spatially chopped intensity on the sky at three frequencies, and in two rows, are recorded. 
The differenced signal is defined as: 

D k = S+ - S- (9) 

where is the signal from each bolometer in the differenced pair. Because of the spacing of 
the photometers, this difference corresponds to a 5' chop on the sky. The subscripts k = 1,2,3 
refer to the frequency bands of 355 GHz (or 273 GHz), 221 GHz and 145 GHz respectively 
in the row that is on the source. The subscripts k = 4,5,6 refer to the same frequency 
set but in the row that is off-source (see Figure 2). In addition to the differenced signal, 
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one bolometer signal from each pair is also recorded, to allow monitoring of common-mode 
signals. These six "single-channel" signals are referred to as Sk, where k is the frequency 
subscript previously described. For example, the difference and single channel at 145 GHz 
of the on-source row will be referenced as D 3 and S3. Both the differences and the single 
channels are sampled at 7 Hz. 

The first step in data analysis is to remove cosmic ray spikes by carrying out a point by 
point differentiation of data from a single scan and looking for large (> 4a) deviations from 
the noise. With a knowledge of the time constant of the bolometer and the height of the 
spike, we can make a conservative estimate of how much data is contaminated and exclude 
that data from our analysis. To account for the effect of the bolometer time-constant, which 
is < 100 ms, we flag a region 100 ms x \n(V m /a) before, and 250 ms x \n(V m /a) after the 
spike's maximum where V m is the height of the spike. The data are then combined into 3 
second bins each containing 21 samples and covering a region equal to 0'75cos<5 on the sky. 
Bins with 11 or more contaminated samples are excluded from further analysis. Less than 
1% of the data are discarded due to cosmic ray contamination. 

4. Analysis of Calibrated Data 

Each cluster data set typically comprises several hundred drift scans, as summarized in 
Table 2. Once the data have been despiked, binned and calibrated, we need to extract the SZ 
signal from the data, and at the same time, obtain an accurate estimate of the uncertainty. 
For most observing conditions atmospheric emission dominates the emission from our source. 
Even under the best conditions the NEF for the 145 GHz channel, which is least affected 
by the atmosphere, is rarely below 100 mJy s 1 / 2 at the signal frequencies of interest, which 
are ~ 10 mHz, while the signals we are trying to measure are ~ 50 mJy. The higher 
frequency channels are progressively worse. In order to improve our signal/noise, we make 
use of our ability to simultaneously measure the sum of the SZ signal and the atmosphere at 
three different frequencies. The different temporal and spectral behavior of the atmosphere, 
compared to the SZ signal of interest, allows us to clean the data and significantly improve 
the sensitivity of our measurements. 

4.1. SZ Model 

A model for the expected spatial distribution of the SZ signal in each scan is obtained 
by convolving a beam map of a planetary calibrator with the modelled opacity of the cluster. 
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Beam shapes are measured by performing raster scans of a planetary calibrator and recording 
the voltage response of the detectors, Vk(9,<f>) where 9 is measured in the direction of RA 
and in the direction of declination. We approximate the electron density of the cluster 
with a spherically-symmetric isothermal (3 model (Cavaliere & Fusco-Femiano 1976, 1978): 



n e (r) = n e0 



l -3/3/2 



1 + 



(10) 



where r is distance from the cluster center, and (3 and r c are parameters of the model. By 
integrating n e along the line of sight, the cluster optical depth: 



t(9,<P)=t 



1 + 



+ 



91 



(l-3/3)/2 



is obtained, where linear distance r has been replaced with angles on the sky, 9 and 0. 
The model parameters of the intra-cluster gas, (3 and 9 C , for each cluster are taken from 
the literature and are listed in Table 4 with associated references. We can now calculate a 
spatial model, m k (9), for each cluster: 



m k {9) 




V k {9',4>')T{9 l -9,<t> 1 ) 



v., 



peak 



T 



d9'd(f)' 



(12) 



that has units of steradians, and is calculated at 0'.05 x cos5 intervals for a given offset, 9, 
in right ascension from the cluster center. We calculate our SZ model by multiplying the 
source model by thermal and kinematic band-averaged spectral factors given by: 



m e c 2 JV(x,T e ) x f k (x)dx 



S fk(x)d. 



and 



K k = -I x 



m e c 
~kT~ e 



x 



1 / h(x,T e ) x f k (x)dx 

X 



(13) 



(14) 



f f k (x)dx 

where the spectral functions ty(x,T e ) and h(x,T e ) were previously defined in §2, and f k (x) 
is the spectral response of channel k. The vector n„ is a unit vector in the direction of the 
cluster peculiar velocity. The quantities T k x m k {9) and K k x m k {9) are then the SZ models 
for the expected responses of frequency band A; to a scan across a cluster of unity central 
comptonization, y , with a radial component to the peculiar velocity, v p , of lkms" 1 . The 
calculated SZ model is then combined into 0'.75 x cos 5 bins to match the binned SuZIE II 
data, so that we define T k x m k (9i) as the thermal SZ model in channel k for the right 
ascension offset 9 of bin number i. 
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4.2. Removal of Residual Atmospheric Signal 

There are two sources of residual atmospheric noise in our data, with different temporal 
spectra. The first is incomplete subtraction of the signal that is common to each beam be- 
cause of the finite common mode rejection ratio (CMRR) of the electronic differencing. This 
effect is minimized by slightly altering the bias, and thus the responsivity, of one of the two 
detectors that form a difference. This trimming process is carried out at the beginning of an 
observing campaign and is usually left unchanged throughout the observations. The second 
is a fundamental limitation introduced by the fact that the two beams being differenced 
pass through slightly different columns of atmosphere; consequently there is a percentage of 
atmospheric emission that cannot be removed by differencing. While both signals originate 
from the atmosphere, their temporal properties are quite different and are accordingly re- 
moved differently in our analysis. In what follows we denote each frequency channel with 
the subscript k, each scan with the subscript j and each bin within a scan with a subscript 
i. In this way the difference and single channel signals at 145 GHz from scan j and bin i are 
D 3ji and S 3ji . 

The residual common mode signal from the atmosphere in the difference channel Dj-ji 
is modelled as proportional to the signal from the corresponding single channel. We define 
our common mode atmospheric template, C^, as Ckji = S^ji- Because the single channels 
contain a small proportion of SZ signal, there is potential to introduce a systematic error by 
removing true SZ signal. However, the effect is estimated at less than 2% (see §5.3.2). 

To model the residual differential signal from the atmosphere we construct a linear 
combination of the three differential channels in a single row which contains no thermal or 
kinematic SZ signal. For the on-source row we define our differential atmospheric template, 

dbS. 

Aji = aD^ + 7-D 2 ji + D Zji (15) 

with a similar definition for the off-source row. The coefficients a and 7 are chosen to 
minimize the residual SZ flux in A^. We describe the construction of this template in detail 
in appendix A and list the values of a and 7 used for the on-source row observation of 
each cluster in Table 5. Removing atmospheric signal in this way significantly increases our 
sensitivity; however it has the disadvantage of introducing a correlation between different 
frequency channels which must be accounted for. Also this model is dependant on the cluster 
parameters used, and is subject to uncertainties in the temperature, and spatial distribution 
of, the cluster gas. We quantify the uncertainty in the final result that this produces in 
section 5.3.2. 

In addition to the atmospheric signals, we also remove a slope, 6^, and a constant, a^, 
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such that our "cleaned" signal is then: 

Xkji = Dkji — dkj — ibkj — (ekj x Ckji) — (fkj x Aji) (16) 
Noting that since we remove a best-fit constant offset this implies £V X fe jj 0. 

4.3. Determination of the Cluster Location in the Scan 

Variations in the location of the cluster center with respect to the nominal pointing 
center defined in Figure 2 can be caused by differences in the location of the X-ray and SZ 
peaks, and by CSO pointing uncertainties. The latter are expected to be less than 10". To 
determine the true cluster location we first co-add all of the scans for a single cluster, as 
described below, then fit the data with the SZ model described above, allowing the source 
position to vary. Note we are only able to constrain the location in right ascension; the 
effects of pointing uncertainties are discussed further in §5.3.3. 

Following Holzapfel et al. (1997a), we define X ki , the coadded signal at each location, 

i, as: 

y _ ggi^WgMgj .... 

where N s is the number of scans, and each scan is weighted according to its root-mean-square 
(RMS) residual defined as: 

RMSL = %^ f (18) 

where iVj, is the number of bins in a single scan. The uncertainty of each bin in the co-added 
scan, is estimated from the dispersion about the mean value weighted by the RMS^ of each 
scan, 



Y:% (x kt - x kJt y/RMS 2 



(19) 



\ (n s - 1) Ef=i i/RMS^ 

This expression provides an unbiased estimate of the uncertainty associated with each bin. 

The on-source row at v ~ 145 GHz (k = 3) provides the highest sensitivity measurement 
of the cluster intensity, and so it alone is used to fix the cluster location. The co-added data 
are fit to a model that includes an offset, a, a slope, b, and the SZ model, where the cluster 
location, RA of j set and the central comptonization, y , are allowed to vary. For each set of 
parameters, we can define x 2 as: 

x 2 = ~ {yo x T 3 x m 3 (9j - RApffset)} - a - ib] ^ 



=1 



'3i 
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To determine the best-fit model all four parameters (a,b,yo, and the RA g- sct ) are allowed 
to vary while the x 2 is minimized. Here we are making the assumption that the measured 
SZ emission in the 145 GHz band is entirely thermal. We are not yet concerned with distin- 
guishing thermal from kinematic SZ emission because at this stage our goal is only to fit the 
location of the cluster. The cluster locations determined in this way are listed in Table 6. 
All of the clusters lie within 30" of the nominal pointing center, and in most cases the cluster 
is located at the pointing center, within our experimental uncertainty. 

For several reasons we use the coadded method only to determine the cluster location, 
not to determine yo and v p . These reasons include one pointed out in Holzapfel et al. (1997a), 
which is that if the source contributes significantly to the variance of each scan, then the 
RMS given by equation (18) will be biased. This does not affect the determination of the 
cluster location. Although we could correct this bias by subtracting the best-fit model from 
the data prior to estimating the RMS, in an iterative fashion, there is another more serious 
complication to the coadded data - that of correlations between the bins in the co-added scan 
produced by the presence of residual atmospheric noise in the data, and by the atmospheric 
removal process itself. Neither the bias of the RMS or the correlations in the coadded scan 
affect the determination of the cluster center, but they do need to be correctly accounted 
for in the determination of the SZ parameters. 



4.4. Individual Scan Fits for Comptonization and Peculiar Velocity 

To generate an unbiased estimate of the SZ parameters we fit for comptonization and 
peculiar velocity in all three frequency channels simultaneously, using the cluster central po- 
sition determined from the coadded data. Following Holzapfel et al. (1997a,b) and Mauskopf 
et al. (2000), we fit the data on a scan- by-scan basis to estimate the uncertainty in the fitted 
parameters, because we expect no scan-to-scan correlation in the noise. While unbiased and 
producing satisfactory results, this method is not formally optimal. An alternative would be 
to calculate, and then invert, the noise covariance matrix for the data set. However, because 
of the high degree of correlation in the raw data, this technique has not been found to yield 
stable solutions. 

We begin again with the de-spiked, binned, calibrated data defined in section 3.5. We 
fit the data vector from each scan with a slope, a constant, the model for residual common- 
mode and differential atmospheric signals and an SZ model with thermal and kinematic 
components. Within each scan we allow the slope, constant, and atmospheric coefficients to 
vary between frequency channels, but we fix the comptonization and peculiar velocity to be 
the same at each frequency. The residual signal left after removal of all modelled sources of 
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signal is then: 

Rkji = Dkji—akj —ibtj — (e k j x C k ji) — (fukj x Aji) — y Q jT k m k (9i— RA offset ) — {yov p ) jK k rrik (6>j — RA offset ) 



where the offset terms, b kj are the slope terms, and e k j (f k j) are the coefficients 

that are proportional to the common-mode (differential-mode) atmospheric signal in fre- 
quency channels k — 1, 2, 3. The SZ-model parameters yoj and (yov p )j are proportional to 
the magnitude of the thermal and kinematic components in each frequency channel. The 
common-mode and differential atmospheric templates, C k ji and A,,, are constructed us- 
ing the method described in section 4.2. The thermal and kinematic SZ model templates, 
T k m k (6i — RA offsct ) and K k m k (6i — RA offsct ), are described in section 4.1. 

The best-fit model of scan j is then determined by minimizing the x 2 , which is defined 

as: 



is the mean squared of the residual signal after removal of the best-fit model. This has 
to be an iterative process because we cannot correctly calculate the best fit model and its 
associated uncertainty until we know the RMS of the residual signal with the best-fit model 
removed. As a first guess we use the RMS of the raw data, and upon each iteration afterwards 
calculate the RMS with the best-fit model removed from the previous minimization. This 
process is continued until the best-fit values for y j and (yoV p )j vary by less than one part 
in a million, a condition usually met by the third iteration. We calculate the uncertainty of 
yoj and (y v p )j, which we define as a y j and (?{ yVv )j, from the curvature of x] m the region of 
the minimum. 



(21) 




(22) 



where 




(23) 



4.5. Likelihood Analysis of Individual Scan Fits 



From the individual scan fits for comptonization and peculiar velocity we next define a 
symmetric 2 by 2 covariance matrix, S, defined by 
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where yoj, (yov p )j, a y j and <7( y v p )j are determined for scan j from the minimization of the x] 
defined in equation (22). The quantities < yo > and < y v p > are weighted averages of the 
individual scan fits for the thermal and kinematic SZ components, and are defined as: 

< yo >= -/ ]J ° 1 ' (27) 
Ef iK 

< vov P >= ^ 2 (28) 

These weighted averages are unbiased estimators of the optical depth and peculiar velocity. 
Having calculated the covariance matrix we define the likelihood function for our model 
parameters v p and y as: 

T 

L(v ) = 1 ex f _ I ( < yo > -yo \ s _i ( < yo > -yo \ 

[v p , y ) ( 2 7t)|S| 1 /2 exp [ 2 V < yo^ P > -yo x v p J \< y v p > -y x v p J 

(29) 

The likelihood is calculated over a large grid in parameter space with a resolution of Av p = 
25kms~ 1 and Ay = 10~ 6 . As an example, the combined likelihood for the measurements of 
MS0451 in November 1996, 1997, and 2000 is shown in Figure 3. The degeneracy between 
a decreasing comptonization and an increasing peculiar velocity is a general characteristic 
of the likelihood function of each cluster. The 1-a uncertainty on each parameter is then 
determined using the standard method of marginalizing the likelihood function over the other 
parameter. The results for each cluster are shown in Table 6. 

A comparison of the numbers in Table 6 for A1835 with earlier results reported in 
Mauskopf et al. (2000) that gave a y-parameter of (6.3 ±1.8) x 10~ 4 and a peculiar velocity 
of 900 ± 1500 kms" 1 shows that the agreement is good. 



4.6. Spectral Plots for Each Cluster 



Figure 4 plots the best-fit SZ spectrum for each cluster with the SuZIE Il-determined 
intensities at each of our three frequencies overlaid. Note, these plots are for display purposes 
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only to verify visually that we do indeed measure an SZ-type spectrum. Although the values 
of the intensity at each frequency are correct, the uncertainties are strongly correlated. 
Consequently these intensity measurements cannot be directly fitted to determine SZ, and 
other, parameters. This is why we use the full scan-by-scan analysis described in the previous 
section. 

In order to calculate the points shown in Figure 4 we calculate a new coadd of the data 
at each frequency after cleaning atmospheric noise from the data. We define the cleaned 
data set, V^.^, as: 

Ykji = D k ji — a k j — ib k j — (e k j x C k ji) — (fkj x A k ji) (30) 

with the best-fit parameters for a k j, b k j, e k j, and f k j determined from equation (22). This 
cleaned data set can now be co-added using the residual RMS defined in equation (23) as a 
weight, such that: 

} kj 



Ef=i*WRMS 



Y* = ^ ,' ' (31) 



Unlike equation (18) used in §4.3, this calculation of the RMS is not biased by any contri- 
bution from the SZ source. The uncertainty of each co-added bin, o~ k i, is determined from 
the dispersion about the mean value, Y ki , weighted by the RMS^ of each scan, 



0~ki 



ggi (Yki - n^VRMs^. 
M (N. ~ 1) Eji i/RMS^ 



The best-fit central intensity, I k , for each frequency band is then found by minimizing 
the xt °f ^ e ^ to the co-added data, where xl 1S defined as follows: 

2 _ [Yki — h x m ki (6j — RA offsct )] 

Xk - T2 y 66 ) 

i=l ki 

We calculate confidence intervals for I k using a maximum likelihood estimator, L(I k ) oc 
exp(— In Figure 5, we show co-added data scans for the November 1997 observations 
of MS0451 for all three on-source frequency bands. The best fit intensity at each frequency, 
and the 1-a error bars are also shown. 

In Figure 6, we show the spectrum of MS0451 measured during each of the three observ- 
ing runs, and the averaged spectrum. This figure, and the best fit parameters (determined 
from the scan-by-scan fitting method) shown in Table 6 indicate that there is good consis- 
tency between data sets taken many months apart. 
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In order to demonstrate the value of our atmospheric subtraction procedure, we have 
repeated our analysis for the MS0451 November 2000 data both with and without atmo- 
spheric subtraction. The derived fluxes from the coadded data are shown in Table 7. The 
improvement in the sensitivity, especially at 220 and 355 GHz, is substantial. 



5. Additional Sources of Uncertainty 

The results given in Table 6 do not include other potential sources of uncertainty in 
the data, such as calibration errors, uncertainties in the X-ray data, and systematic effects 
associated with our data acquisition and analysis techniques. We now show that these un- 
certainties and systematics are negligible compared to the statistical uncertainty associated 
with our SZ measurements. Astrophysical confusion is considered separately in §6. 



5.1. Calibration Uncertainty 

To include the calibration uncertainty, we use a variant of the method described in 
Ganga et al. (1997). A flux calibration error can be accounted for by defining a variable, 
Gfc, such that the correctly calibrated data is D k j i = G k x D k ji- We further assume that 
the calibration error can be broken down into the product of an absolute uncertainty that 
is common to all frequency bands, and a relative uncertainty that differs between frequency 
bands. In this way we define G k = G abs x G r k el with the assumption that both G abs and 
G r k el can be well-described by Gaussian distributions that are centered on a value of 1. The 
likelihood, marginalized over both calibration uncertainties, is then 

poo poo 3 

L(y , v p ) = / dG abs P(G abs ) / L(y , v p , G abs , G[ el , G r 2 el , G r 3 el ) TT dG r k el P(G r k d ) (34) 
Jo Jo k=1 

We evaluate these integrals by performing a 3 point Gauss- Hermite integration (see Press et 
al. 1992, for example) using the likelihoods calculated at the most-likely values, and the 1-a 
confidence intervals, for G abs and G r k el , which we will now discuss. While the main source 
of absolute calibration error in our data is the ±6% uncertainty in the RJ temperature 
of Mars and Saturn (see Section 3.4), it is not straightforward to label other calibration 
uncertainties as either absolute or relative. Instead we calculate equation (34) assuming 
two different calibration scenarios: one where our calibration uncertainty is entirely absolute 
such that G abs = [0.9, 1.00, 1.1] and G' k l = [1.0, 1.0, 1.0] for all k, and the other which has 
a equal combination of the two with G abs = [0.93, 1.00, 1.07] and G' k l = [0.93, 1.00, 1.07] for 
all values of k. This allows us to assess whether the assignment of the error is important. 
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We have recalculated the best fit y and v p using the MS0451 data taken in November 
2000. We choose this data set because it has some of the lowest uncertainties of any of 
our data sets and consequently we would expect it to be the most susceptible to calibration 
uncertainties. Ignoring the calibration uncertainty, we calculate yo = 3.06!q s! x and 
v p = —300^250 km s -1 from marginalizing the likelihood as described in §4.5. We then let 
Qabs an( j Qrei var y over their allowed range to calculate L(y , v p , G abs , G r k el ) with a param- 
eter space resolution of Ay = 10 -6 and Av p = 25kms _1 in our two different calibration 
scenarios. Assuming only absolute calibration uncertainty we marginalize this likelihood 
over G abs = [0.9,1.00,1.1] and G r k el = [1.0,1.0,1.0] and find that y = 3.02±°j 4 x 10~ 4 
and v p = — 300l{25o km s -1 , values which are completely unchanged from the best fit 
values assuming no calibration uncertainty. Assuming a combination of absolute and rel- 
ative calibration uncertainty we marginalize this likelihood over G abs = [0.93, 1.00, 1.07] 
and G r k el = [0.93,1.00,1.07]. We find new best fit values of y = 3.02±°;|[ x 10" 4 and 
v p = —275^250 km s _1 , again virtually identical to the values obtained assuming no calibra- 
tion uncertainty. Therefore we conclude that for all our clusters the error introduced from 
calibration uncertainty, regardless of source, is negligible compared to the statistical error of 
the measurement. The effects of calibration uncertainties are summarized in Table 8. 



5.2. Gas Density and Temperature Model Uncertainties 

We now account for the effect of uncertainties in the (3 model parameters for the intra- 
cluster gas by fitting our SZ data with the allowed range of gas models based on the \-a 
uncertainties quoted for (3 and 9 C in Table 4. Ideally one would fit the X-ray and SZ data 
simultaneously to determine the best-fit gas model parameters. For several of our clusters 
this has been done using 30 GHz SZ maps by Reese et al. (2002), and we use the values for 
the f3 model derived in this way. SuZIE II lacks sufficient spatial resolution to significantly 
improve on constraints from X-ray data, and so for clusters that are not in the Reese et al. 
(2002) sample, we use the uncertainties derived from X-ray measurements alone. 

Using a similar method to the calibration error analysis in the previous section, we 
assume that the range of allowable gas models can be well-approximated by a Gaussian 
distribution centered around the most-likely value and marginalize the resulting likelihood 
integrals over (3 and core individually using 3 point Gauss-Hermite integration. In reality, 
the gas model parameters f3 and Q core are degenerate and their joint probability distribution 
is not well-approximated by two independent Gaussians. However, this crude assumption 
allows us to show below that this source of error is relatively negligible compared to the 
statistical error of our results. 
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To estimate the effects of density model uncertainties in our sample we study the effect 
on MS0451 because it has one of the least well constrained density models from our sample. 
We find that when the allowable range of uncertainty on (3 and 6 C is included, the best fit 
SZ parameters are yo = 3.04+Qg 3 x 10~ 4 and v p = — 300^1250 km s_1 ' unchanged from the 
values in table 6. Therefore we conclude that the error from density model uncertainties is 
negligible compared to the statistical error of the measurement. 

To estimate the effects of temperature model uncertainties in our sample we again use 
MS0451 because it has one of the least constrained electron temperatures from our sample. 
We again assume that the range of allowable temperatures is well approximated by a Gaus- 
sian distribution centered around the most-likely value and marginalize the resulting likeli- 
hood integrals over T e using 3 point Gauss-Hermite integration. We find yo — 3.05^g ^ xlO" 4 
and v p = — 275^250 km s_1 > virtually unchanged from the best fit values that assume no 
temperature uncertainty. Therefore we conclude that the error from temperature uncertain- 
ties is negligible compared to the statistical error of the measurement. The effects of gas 
model uncertainties are summarized in Table 8. 



5.3. Systematic Uncertainties 

We now consider effects that could cause systematic errors in our estimates of y and 
v p . These include instrumental baseline drifts that could mimic an SZ source in our drifts 
scan, and systematics introduced by our atmospheric subtraction technique. 

5.3.1. Baseline Drifts 

Previous observations using SuZIE II, and the single-frequency SuZIE I receiver, have 
found no significant instrumental baseline (Mauskopf et al. 2000; Holzapfel et al. 1997a,b). 
Baseline checks are performed using observations in patches of sky free of known sources or 
clusters. For the data presented in this paper we also use measurements with SuZIE II on re- 
gions of blank sky. In February 1998 we observed a region of blank sky at 07h40m0 s ;+9 o 30'0" 
(J2000) for a total of ~ 18 hours of integration in exceptional weather conditions. The sky 
strip was 60' in length and was observed in exactly the same way as the cluster observations 
presented in this paper. This data represents the most sensitive measurements ever made 
with SuZIE II and consequently should be very sensitive to any residual baseline signal (the 
data itself will be the subject of a separate paper). We have repeated exactly the analysis 
procedure used to analyze our cluster data with one exception, that we restrict the test 
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source position in the blank sky field to be within O'.l of the pointing center. The source 
positions derived from the cluster data are consistently < 30" from the pointing center indi- 
cating that there is no significant off-center instrumental baseline signal. For the fit, we use a 
generic SZ source model with (3 = 2/3 and 9 C = 20" and find the best-fit flux to the co- added 
data in each on-source frequency band. In our 145 GHz channel we find a best fit flux of 
AI = -1.05 ±2.12 mJy in our 221 GHz channel we find a best fit flux of AI = -3.56 ±3.79 
mJy, and in our 355 GHz channel we find a best-fit flux of AI = —1.91 ± 6.94 mJy. Using 
the November 2000 data from MS0451 as an example, if the blank sky flux measured at 145 
GHz was purely thermal SZ in origin this would correspond to a central comptonization of 
y = 0.10 ± 0.20 x 10~ 4 , while the blank sky flux measured at 221 GHz, assuming r = 0.015, 
corresponds to a peculiar velocity of v p = 635 ± 676 km s" 1 . Therefore we conclude that 
there is no significant systematic due to baseline drifts in any of our three spectral bands. 



5.3.2. Systematics Introduced by Atmospheric Subtraction 

The model that was fitted to each data set, D k ji, as defined in equation (21), included 
common-mode atmospheric signal, Ckji, that was defined to be proportional to the average of 
our single channel signals, S^i- While the single channel signal is dominated by atmospheric 
emission variations, it will also include some of the SZ signal we are trying to detect. This 
can potentially cause us to underestimate the SZ signal in our beam because part of it will 
be correlated with the single channel template. We estimate this effect from the correlation 
coefficients e k j calculated during the minimization of in equation (22). We estimate that 
the total SZ signal subtracted out from our common-mode atmospheric removal is ~ 2%, at 
a level that is negligible compared to the statistical error of our results. 

The construction of a differential atmospheric template can potentially introduce resid- 
ual SZ signal through our atmospheric subtraction routine. We discuss our method to 
construct a differential atmospheric template in appendix A and follow the notation defined 
therein. Residual SZ signal in this template can be introduced through the simplifying as- 
sumption that the a and 7 used in its construction are spatially independent. In addition, 
uncertainties in the electron temperature and density model of the cluster affect how accu- 
rately a and 7 are defined. Below we examine the effects of these two sources of uncertainty 
for the November 2000 observations of MS0451. 

To model the effect of a residual SZ signal in our atmospheric template, A,-j, we re- 
define it by subtracting out the expected residual thermal and kinematic signals, and 
Zj^ (defined in appendix A), binned to match the data set. To calculate Z^ and Zj^ we 
use the values of a and 7 given in Table 5 and assume the comptonization and peculiar 
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velocity values given in Table 6. For all the clusters in our set we find \Z^(9)\ < 5.7 mJy and 
\Zj?(9) | < 1.7 mJy across a scan of the cluster. Using the re-defined atmospheric template we 
then repeat the analysis of the data set and recalculate comptonization and peculiar velocity. 
Using the November 2000 data of MS0451 as an example, we calculate yo = 3.06lo'g3 x 10~ 4 
and v p = —300^1250 km s_1 usm g the method described in section 4.5. Using these values 
for comptonization and peculiar velocity we re-define our atmospheric template as described 
above. We then repeat our analysis routine exactly and calculate y = 3. 021^82 x 10~ 4 an d 
v p = -2501^ km s" 1 . 

The accuracy of the construction of our differential atmospheric template, parameterized 
by the variables a and 7, is limited by our knowledge of each cluster's density model and 
electron temperature. We have calculated a and 7 for each cluster using the best-fit spherical 
beta model parameters 9 C ), and electron temperature (T e ). We recalculate a and 7 using 
the la range of T e ,(3 and 9 C . For the November 2000 observations of MS0451, variations in 
the model parameters of the cluster cause < 1% changes in a and 7. Using the most extreme 
cases of a and 7 we find changes of ±0.01 x 10~ 4 in yo and ±25 km s _1 in v p . Adding the two 
sources of error of differential atmospheric subtraction, discussed in the above paragraphs, 
in quadrature we find an overall uncertainty of 1q 04 x in y and +25 km s" 1 in v p . We 
therefore conclude that uncertainty from differential atmospheric subtraction adds negligible 
error compared to the statistical error of our results. 



5.3.3. Position Offset 

In section 4.3 we allow the position of our SZ model to vary in right ascension and 
determine confidence intervals for this positional offset. However, we do not have the nec- 
essary spatial coverage to constrain our clusters' position in declination, 5. If a cluster's 
position was offset from our pointing center in declination, we would expect the measured 
peak comptonization to be underestimated from the true value. From observations of several 
calibration sources over different nights, we estimate the uncertainty in pointing SuZIE to 
be < 15". The cluster positions that we use are determined from ROSAT astrometry, which 
is typically uncertain by ~ 10 — 15". Adding these uncertainties in quadrature we assign an 
overall pointing uncertainty of AS ~ 20". 

To estimate the effects of pointing uncertainty in our sample we study the effect on 
observations of MS0451 in November 2000. Using the method described in section 4.5, 
which assumed no pointing offset, we calculated yo = 3.06l 'g3 x 10~ 4 and v p = — 300^250 
km s" 1 . We re-calculate the SZ model of MS0451 with a declination offset of 20" from our 
pointing center. Using this SZ model we repeat our analysis routine exactly and calculate 
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yo = 3.18lo|g x 10 -4 and v p = -3001^50 km s_1 - This corresponds to a ~ 4% underestimate 
of the peak comptonization, however there is no effect on peculiar velocity. 

6. The Effects of Astrophysical Confusion 

6.1. Primary Anisotropies 

Measurements of the kinematic SZ effect are ultimately limited by confusion from pri- 
mary CMB anisotropies which are spectrally identical to the kinematic effect in the non- 
relativistic limit. LaRoque et al. (2002) have estimated the level of CMB contamination in 
the SuZlE 11 bands for a conventional Q — 1 (ACDM) cosmology using the SuZIE II beam 
size, at \5y \ < 0.05 x I0~ 4 and \5v pec \ < 380 km s -1 . At present this is negligible compared 
to our statistical uncertainty. 

6.2. Sub-millimeter Galaxies 

Sub-millimeter galaxies are a potential source of confusion, especially in our higher 
frequency channels. All of our clusters have been observed with SCUBA at 450 and 850 /zm 
and sources detected towards all of them at 850/im (Smail et al. 2002; Chapman et al. 2002). 
Because of the extended nature of some of these sources, it is difficult to discern which are 
true point sources and which are, in fact, residual SZ emission (Mauskopf et al. 2000). This 
is especially true when the source is only detected at 850/im. We assume a worst-case - that 
all of the emission is from point sources - and examine the effects of confusion in MS045I, 
AI835, and A2261. We select these clusters because the sources in MS0451 and A2261 have 
SCUBA fluxes typical of all of the clusters in our set, while A1835 has the largest integrated 
point source flux, as measured by SCUBA, of all our clusters. We consider only sources with 
declinations that are within V in declination of our pointing center since these sources will 
have the greatest effect on our measurements. The point sources that meet this criterion 
are shown in Table 9. To model a source observation, we use SCUBA measurements to set 
the expected flux at 850yum, and assume a spectral index, a, of 2 or 3 to extrapolate to 
our frequency bands, where the flux in any band is S v oc u a . In the case of Af835, where 
450 fim fluxes have also been measured, the spectral indices of the detected sources range 
from f .3-2, with large uncertainties. For this cluster we also examine the effect of an index 
of 1.7. 

Note that because our beam is large we cannot simply mask out the SCUBA sources 
from our scans without removing an unacceptably large quantity of data. Instead, for each 
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cluster, the sources are convolved with the SuZIE II beam-map to create a model observation 
at each frequency. This model is then subtracted from each scan of the raw data, and 
the entire data set is re-analyzed. The results are summarized in Table 10. The overall 
effect of sub-millimeter point sources is to increase the measured flux at each frequency by 
10 — 50% of the point source flux at that frequency. The reason that the flux error is less 
than the true point source flux is that spectrally the point sources are not too different 
from atmospheric emission, which also rises strongly with frequency, and so the sources 
are partially removed during our atmospheric subtraction procedure. Because the residual 
point source flux is the same sign at all frequencies, it is spectrally most similar to the 
kinematic effect. Consequently the effect on the final results is to slightly over-estimate 
the comptonization, with the peculiar velocity biased to a more negative value by several 
hundred kms -1 . While the uncertainty this introduces is currently small compared to the 
statistical uncertainty of our measurements, it is a systematic that will present problems for 
future more sensitive measurements and we expect that higher resolution observations than 
SCUBA will be needed in order to accurately distinguish point sources from SZ emission. In 
addition the SCUBA maps cover only 2(3, and so all of the sources that we consider here cause 
a systematic velocity towards the observer. Because of the differencing and scan strategy 
that we use, sources that lie outside the SCUBA field of view can cause an apparent peculiar 
velocity away from the observer which is not quantified in this analysis. We conservatively 
estimate this contribution to be equal in magnitude but opposite in sign to the effect of 
sources within the field of view. In reality, the effects will likely cancel to some degree, 
reducing the overall uncertainty associated with submm sources. 



6.3. Unknown Sources of Systematics 

Finally, in order to check for other systematics in our data, we calculate the average 
peculiar velocity of the entire set of SuZIE II clusters, taking into account the likelihood 
function of each measurement based on the statistical uncertainties only, and find that the 
average is ISOlggQ km/s. This indicates that, at present, our experimental uncertainties 
exceed any systematic errors in our data. 



6.4. Summary 

Table 11 summarizes the effect of all of the known sources of uncertainty in our mea- 
surement of the peculiar velocity of MS0451. We expect similar uncertainties for the other 
clusters in our sample. Other than the statistical uncertainty of the measurement itself, the 
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dominant contribution is from CMB fluctuations and point sources. Note that we do not 
include the calculation of the baseline presented in §5.3.1 because our measurements show 
no baseline at the limit set by astrophysical confusion. 

7. Checking for Convergence of the Local Dipole Flow 

Table 12 summarizes the current sample of SuZIE II measurements of peculiar velocities 
and includes two previous measurements made with SuZIE I of A2163 and A1689 (Holzapfel 
et al. 1997b). Note that in each case we have assumed the statistical uncertainty associated 
with each measurement, then added in quadrature an extra uncertainty of ±750 km s -1 to 
account for the effects of astrophysical confusion from the CMB and submm point sources, 
based on the estimates derived in §6. These confusion estimates are larger than presented in 
Holzapfel et al. (1997b), mainly because of the more recent data on submm sources. Conse- 
quently we use values for the peculiar velocities of A2163 and A1689 that have uncertainties 
that are somewhat higher than those previously published. 

The locations of the clusters on the sky are shown in Figure 7. The figure also shows 
the precision on the radial component of each cluster peculiar velocity, plotted against the 
redshift of the cluster. The cross-hatched region shows the region of redshift space that 
has been probed by existing optical surveys (see below). We now use this sample of cluster 
peculiar velocities to set limits to the dipole flow at these redshifts. The clusters in the 
SuZIE II sample all lie at distances > 500/1" 1 Mpc, where the flow is expected to be < 50 
km s^ 1 (see below). Of this sample, we have taken the 6 that lie in the range z = 0.18-0.29 
and used them to set limits on the dipole flow of structure at this redshift. This calculation 
gives a useful indication of the abilities of SZ measurements to probe large-scales motions 
and is also the first time that such a measurement has been made using SZ results. 

Flows that are coherent over large regions of space probe the longest wavelength modes 
of the gravitational potential and provide a test of matter formation and evolution in the 
linear region. The average bulk velocity, Vb, of a region of radius, R, is predicted to be: 

rr2ol.2 poo 

( V H R )) = J P(k)W 2 (kR)dk (35) 

(see Strauss & Willick 1995, for example), where Vt m is the present-day matter density, H is 
the Hubble constant, P(k) is the matter power spectrum and W(kR) is the Fourier transform 
of the window function of the cluster sample. The shape of the window functions depends 
on details of the cluster sample such as whether all of the sphere is sampled, or whether the 
clusters lie in a redshift shell (Willick 1999). In all models the flow is expected to converge 
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as a function of increasing R and at z > 0.18 it should be less than 50kms _1 . Of course, 
equation (35) is strictly valid only in the local universe. At higher redshifts, the rate of 
growth of fluctuations must be accounted for, causing the expected bulk flow to decrease 
even more quickly. 

All measurements of bulk flows to date have used peculiar velocities of galaxies, deter- 
mined by measuring the galaxy redshift and comparing it to the distance determined using a 
distance indicator based on galaxy luminosities, rotational velocities, or Type la supernovae. 
These methods have yielded bulk flows consistent with theory out to a distance of 60/i _1 Mpc 
(see Courteau & Dekel 2001, for reviews of the experimental situation). At larger distances, 
the situation is less clear. Some null measurements do seem to confirm that the flow con- 
verges (Colless et al. 2001; Dale & Giovanelli 2000), but there are a number of significantly 
larger flow measurements that have not yet been refuted and that are quite discrepant with 
one another in direction (Willick 1999; Hudson et al. 1999; Lauer & Postman 1994). The 
directions of these flows are shown in Table 13. 

We define a bulk flow as V B = (V B , l B , b B ) where V B is the velocity of the bulk flow, l B 
is the galactic longitude and b B is the galactic latitude of the flow direction. We can then 
calculate the likelihood of these parameters as: 

^(V B ) = n^( V B-^) ( 36 ) 

i 

where Vb-^ is the component of the bulk flow in the direction of the ith cluster, and the 
likelihood of the flow, given the data, Li, is determined from the SuZIE II data. In order to 
account for the effects of astrophysical confusion, we convolve the likelihood function for the 
peculiar velocity of each cluster with a gaussian probability with a = ±750kms _1 before 
calculating the bulk-flow likelihood. As expected, we do not detect any bulk flow in our 
data. Figure 8 shows the 95% confidence limit to Vb as a function of location on the sky. 

Because our clusters sparsely sample the peculiar velocity field, our upper limits are 
tighter in some directions than in others. For example, we have also used our data to set 
limits on bulk flows in the direction of the CMB dipole which is taken to have coordinates 
(l,b) = (276° ± 3°, 33° ± 3°) (Kogut et al. 1993). We find that at 95% confidence the flow 
in this direction is < 1410 km s" 1 . The limits towards other directions for which optical 
measurements have yielded a flow detection are shown in Table 13. 
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8. Conclusions 

We have used millimetric measurements of the SZ effect to set limits to the peculiar 
velocities of six galaxy clusters. By making measurements at three widely separated frequen- 
cies we have been able to separate the thermal and kinematic SZ spectra. Moreover, because 
our measurements in these bands are simultaneous, we have been able to discriminate and 
remove fluctuations in atmospheric emission which dominate the noise at these wavelengths. 
These observations have allowed us to make the first SZ-determined limits on bulk flows. 
In certain directions the limits are approaching the level of sensitivity achieved with optical 
surveys at much lower redshifts. Because our sky coverage is not uniform, our sensitivity to 
bulk flows varies greatly over the sky. We are continuing the SuZIE II observation program 
with the aim of obtaining a larger, more uniform cluster sample to improve our bulk flow 
limits. 

The precision with which SuZIE II can measure peculiar velocities is limited not only 
by the small number of detectors, but also by atmospheric and instrumental noise, and 
astrophysical confusion from submillimeter galaxies. In the last few years the potential of 
SZ astronomy has been realized and new telescopes equipped with bolometer arrays that will 
have hundreds to thousands of pixels are planned (Carlstrom, Holder & Reese 2002; Staggs 
& Church 2001). Our measurements demonstrate the feasibility of using the SZ effect to 
measure cluster peculiar velocities, but highlight some issues that need further investigation: 

1. The atmosphere will continue to be an issue for ground-based measurements. We have 
demonstrated that wide, simultaneous spectral coverage can significantly diminish, but 
not completely remove, this source of noise. 

2. Confusion from submillimeter point sources can lead to a systematic peculiar velocity 
measurement. Further work is needed to determine the best strategy for identifying 
and removing this contaminant, especially for experiments that will map a large area 
of sky with lower angular resolution than SuZIE II, such as the High Frequency In- 
strument (HFI) that will operate at frequencies of 100-850 GHz as part of the Planck 
satellite (Lamarre et al. 2000). Although high-resolution measurements of clusters 
at millimeter wavelengths will be possible with future instruments such as ALMA, it 
may not be feasible to use this method to subtract astrophysical contamination from 
a sample of the size needed (of order 1000+ velocities) to effectively probe large-scale 
structure with SZ-determined peculiar velocities. Instruments such as SuZIE II will be 
invaluable for investigating removal techniques based on spectral rather than spatial 
information, especially if more frequency bands are incorporated into the instrument. 

The SuZIE program is supported by a National Science Foundation grant AST-9970797 
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National Science Foundation Graduate Research Fellowship and a Stanford Graduate Fel- 
lowship. This work was partially carried out at the Infrared Processing and Analysis Center 
and the Jet Propulsion Laboratory of the California Institute of Technology, under a con- 
tract with the National Aeronautics and Space Administration. We would like to thank the 
CSO staff for their assistance with SuZIE observations. The CSO is operated by the Caltech 
Submillimeter Observatory under contract from the National Science Foundation. 



A. Creating an Atmospheric Template 

SuZIE makes simultaneous multi-frequency observations of a source. Besides the obvious 
advantage of instantaneously measuring the spectrum of the SZ source, this method permits 
atmospheric noise removal through spectral discrimination of the atmospheric signal. We 
realize this by creating an atmospheric template with no residual SZ signal for each row of 
photometers by forming a linear combination of the three differential frequency channels in 
that row on a scan by scan basis. We define this atmospheric template, A,-(0), as: 

Aj(9) = aD l3 {6) + rfD v (0) + D 3j (9) (Al) 

where the coefficients a and 7, are chosen to minimize the residual SZ signal in Aj(0). If we 
define F k sz (9) as the SZ signal in channel D k {9) this implies that the residual SZ signal in 
our atmospheric template, which we define as Z k z (9), is then 

Z s k z {9) = aF 1 sz (9) + ^Fi z {9) + F 3 SZ (9) (A2) 

where the SZ flux in each channel includes contributions for both the thermal and kinematic 
effects, such that, F k z {9) = F k Z)T {9) + F k z ' K (9). The thermal SZ flux in each channel, 
using the notation of §4.1, is more precisely: 

F k Z ' T (9) =y x T k (T e ) x m k {9) (A3) 

with yo being the central comptonization of the cluster, T k (T e ) being defined in equation (13), 
and m k {9) defined in equation (12). Similarly, the kinematic SZ flux in each channel is 

F* Z ' K {0) = Vov P x K k (T e ) x m k {9) (A4) 

with v p being the peculiar velocity of the cluster, and K k (T e ) defined in equation (14). 

Using this notation we rewrite our expression for the residual SZ signal in the atmo- 
spheric template as 

Z S k Z {6) = Z T k {9) + Z«{9) (A5) 
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with: 



z T M 



y [a71(T e )mi(0) + 7 T 2 (T e )m 2 (fl) + T 3 (T e )m 3 (9)] 
-y v p [aK^m^O) + 1 K 2 (T e )m 2 (6) + K 3 (T e )m 3 (e)} 



(A6) 
(A7) 



Solving for Zj? z (9) = requires a position-dependent solution for a and 7 because rrik(9) is 
not constant between frequency bands. For simplicity, we use the peak values of rrik{6) to 
calculate a and 7 for Zf z = 0. The error introduced from this assumption is discussed in 
Section 5.3.2. As an aside, this template should also be free of primary CMB anisotropy, on 
spatial scales similar to our clusters, because the kinematic SZ effect is spectrally indistin- 
guishable from a primary CMB anisotropy. 
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Table 1. SuZIE Frequencies and Beam Sizes 



1996 1997-present 

fwhm Ti V Au fwhm n 

Channel 3, [GHz] [GHz] arcmin arcmin 2 [GHz] [GHz] arcmin arcmin 2 



270/350 GHz 



°10 


274.0 


32.6 


1 


.35 


1 


.85 


355.1 


30.3 


1 


.50 


2.17 




272.8 


28.8 


1 


.30 


1 


.79 


354.2 


31.1 


1 


.54 


2.42 


st 


272.9 


30.3 


1 


.25 


1 


.65 


356.7 


30.8 


1 


.57 


2.24 


Sn 


272.3 


30.6 


1 


.30 


1 


.76 


354.5 


31.5 


1 


.60 


2.58 


221 GHz 
























&20 


222.2 


22.4 


1 


.30 


1 


.64 


221.4 


21.8 


1 


.36 


1.90 


^20 


220.0 


23.1 


1 


.25 


1 


.67 


220.5 


23.8 


1 


.40 


2.08 


^21 


220.7 


24.8 


1 


.15 


1 


.53 


221.7 


22.9 


1 


.47 


2.02 


^21 


220.0 


24.1 


1 


.20 


1 


.54 


220.8 


21.7 


1 


.40 


2.18 


145 GHz 
























^30 


146.1 


19.9 


1 


.45 


2 


.14 


145.2 


18.1 


1 


.50 


2.35 


^30 


146.4 


21.3 


1 


.40 


2 


.16 


144.6 


18.6 


1 


.64 


2.73 


S31 


145.0 


20.7 


1 


.30 


1 


.86 


145.5 


18.2 


1 


.60 


2.51 


^31 


144.4 


19.5 


1 


.35 


1 


.91 


144.9 


17.4 


1 


.60 


2.73 



a The notation is of the form Sf reqirow and the ± sign refers to the sign of the channel in 
the differenced data. 
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Table 2. Summary of SuZIE observations 







R.A. a 


Decl. a 




Total 


Accepted 


Integration 




Source 


z 


(J2000) 


(J2000) 


Date 


Scans 


Scans 


Time (hours) 


Ref. 


A 99R1 


u.zz 


1 7 99 97 f\ 


_l_Q9 n7 VI 1 


A/Tar QQ 


1 AA 


loo 


A A 


1 
1 


A2390 


0.23 


21 53 36.7 


+ 17 41 43.7 


Nov 00 


147 


128 


4.3 


1 


ZW3146 


0.29 


10 23 38.8 


+04 11 20.4 


Nov 00 


151 


131 


4.4 


1 


A1835 


0.25 


14 01 02.2 


+02 52 43.0 


Apr 96 


638 


577 


19.2 


1 


C10016+16 


0.55 


00 21 08.5 


+ 16 43 02.4 


Nov 96 


304 


277 


9.2 


2 


MS0451. 6-0305 


0.55 


04 54 10.8 


-03 00 56.8 


Nov 96 


405 


348 


11.6 


2 


// 








Nov 97 


236 


211 


7.0 




// 








Nov 00 


375 


306 


10.2 




MS0451. 6-0305 








Total 


1016 


865 


28.8 





References. — (1) Ebeling et al. (1998); (2) Gioia & Luppino (1994) 

a Units of RA are hours, minutes and seconds and units of declination are degrees, arcminutes and 
arcseconds 



Table 3. Break down of the calibration uncertainties 



Source 


Uncertainty (%) 


Detector non-linearities 


6 


Planetary temperature 


6 


Atmospheric r 


2 


Spectral response 


1 


Beam uncertainties 


5 


Total 


10 
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Table 4. IC gas temperatures beta model parameters 





kT e 








Cluster 


(keV) 




(arcsec) 


Ref. 



A2261 


o OO+0.37 
8 - 82 -0.32 


0.516^ 




1;2;2 


A2390 


iu.io_ 099 


0.67185 


52.018:8 


1;3;3 


Zw3146 


fi 41 +a26 


u - ' ^-0.0 


13.018:8 


1;3;3 


A1835 




u.jyj_ 005 


i 9 9+O.6 
iz - z -0.5 


1;2;2 


C10016 


7 crcr+0.72 


o.749i8;8^ 


42.3±|:S 


4;2;2 


MS0451 


10.4l i 


o.806±8;gti 


34.711:1 


5;2;2 



References. — (1) Allen & Fabian (1998), (2) Reese 
et al. (2002), (3) Ettori & Fabian (1999), (4) Hughes 
& Birkinshaw (1998), (5) Donahue (1996) 



Table 5. Differential Atmospheric Template Factors 



Cluster 


a 


7 


A2261 a 


0.6366 


-1.4490 


A2390 a 


0.6563 


-1.4721 


Zw3146 a 


0.6428 


-1.3808 


MS0451 a (Nov97) 


0.6770 


-1.3933 


MS0451 a (Nov00) 


0.6449 


-1.4286 


Al835 b 


1.2133 


-2.1745 


C10016 b 


1.2353 


-2.2331 


MS0451 b (Nov96) 


1.2246 


-2.2076 



a High frequency channel was 355 
GHz 

b High frequency channel was 273 
GHz 
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Table 6. Summary of Results 



Cluster 


Date 


A RA (arcsec) 


S/o x 10 4 


v pec (km s x ) 


A2261 

A2390 

Zw3146 

A1835 

C10016 


Mar99 
NovOO 
NovOO 
Apr96 
Nov96 


i q+11.0 
1 -°-11.0 
o p;+14.3 
z - d -14.3 
io n+23.9 
io.y_ 22 .9 

99 c+12.0 
_9 q+25.0 


7-58±I;|| 
1.671^ 
3.601^ 
7.541^ 
2.98l£jg 


-140011™ 
+195012$! 

oou -1875 
99^+1650 
— 1225 

-40501??°° 


MS0451 


Nov96 
Nov97 
NovOO 
Combined Fit 


1 cr+16.0 
1 - J -16.0 

22.011° 
-25 +17 '° 


9 i c+0.72 

3-06l°;l 
2.80l°J 2 


_i_oc;n+ 6625 
+ oou -2925 

+ 16501^1 
"300111° 
+75011?°° 



Table 7. The Effects of Atmospheric Subtraction on Derived Fluxes for MS0451 Nov 2000 

Data 



Frequency Flux (MJy sr x ) 

(GHz) With Atm. Subtraction Without Atm. Subtraction 



355 
221 
145 
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Table 8. Effects of Calibration and IC Gas Model Uncertainties on MS0451 (Nov 2000) 



Uncertainty 


yo x io 4 


v p (km 


s- 1 ) 


Statistical Uncertainty 


°- uu -0.83 


-300^ 


-1950 
-1250 


Calibration a (Absolute Only) 


o nO+0-84 
°- uz -0.83 


_OQQ + 1950 

,3UU -1250 


Calibration a ( Equal Absolute and Relative) 


q piQ+0.83 
°- uz -0.84 


-275^ 


-1950 
-1250 


IC Density Model 


o Q4+0.84 
°- u ^-0.83 


-3001^ 


IC Gas Temperature 


o qp:+0.82 
°- ud -0.82 


97 r+1950 
Z ' J -1250 



a See text for details 



Table 9. SCUBA sources towards MS0451, A1835 and A2261. 



Source Coordinates Flux (mJy) 



Cluster 


Source 


RA a (J2000) 


Dec a (J2000) 


850/iin 


450/zm 


Ref. 


A2261 


SMMJ17223 + 3207 


17 22 20.8 


+32 07 04 


17.6 + 3.9 




1 


A1835 


SMMJ14009 + 0252 


14 00 57.7 


+02 52 50 


14.5 + 1.7 


33 + 7 


2 




SMMJ14011 + 0252 


14 01 05.0 


+02 52 25 


12.3 + 1.7 


42 + 7 


2 




SMMJ14010 + 0252 


14 01 02.3 


+02 52 40 


5.4 + 1.7 


20 + 7 


2 


MS0451 


SMMJ04542 - 0301 


04 54 12.5 


-03 01 04 


19.1 + 4.2 




1 



References. — (1) Chapman et al. (2002); (2) Smail et al. (2002) 

a Units of RA are hours, minutes and seconds and units of declination are degrees, arcmin- 
utes and arcseconds 
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Table 10. MS0451, A2261, and A1835 Sub-mm Galaxy Confusion 





Spectral 




Flux/mJys 








Cluster 


Model 


145GHz 


221GHz 


273GHz 


355GHz 


Vo x 10 4 


(km s _1 ) 


MS0451 


No Source 3, 
a = 3 
a = 2 


-24.2^1 
-24.9111 
— 26.3±t;i 


_^ 4+6-7 

-7 ^+ 6 - 7 
'•°-6.7 

-9-51S:? 




+55.41H? 

, n 7 +15.7 
-l-Ol./ _ 15 8 

+49.4111?;? 


o n6 +0 - 83 

O.UO_ 083 

2 Q2 +0 - 81 

z - yz -0.83 

2.871°; 8 3 


qqq+1950 
, 9 p-+2150 

+3501^ 


A2261 


No Source 3, 
a = 3 
a = 2 


-44.5^1 
-45.9+t! 


+n.oill:i 
+8.9iltl 
+7.9lltl 




+113.9l^; 8 7 
+ 110.0™ 
+108.911; 8 


7.58l|:|* 
7.341H? 


-140011™ 
-1250ll 8 ^ 
-1150ll 8 t° 


A1835 


No Source 3, 
a = 3 
a = 2 
a = 1.7 


— 36.3±|;f 
-37.011;? 
-38.01^ 
-38.3l|;? 


+2.9ll;l 

+0.9l|« 
-0 ^ +8 - 6 

u -°-8.6 

-1 o +8 - 6 

i - u -8.6 


+38.2+H? 
+35.21H? 
+33.6ll| 6 7 
+33.011^ 




7 (J4+1-60 

' - d ^-1.61 
7 on+1.61 

7 9C+I.6I 

7 oc+1.61 


-200111° 

, A+1750 

~r u -1300 

, lnn +1800 
-t-iuu_ 1300 

,i c-q+1800 

-|-±ou_ 1300 



This row assumes that there are no sources in the data 
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Table 11. Comptonization and Peculiar Velocity Uncertainties for MS0451 (Nov 2000) 



Uncertainty 


Ho x 10 4 


v p (km s x ) 


Statistical: 


o n fi +0.83 
°- uu -0.83 


QQQ + 1950 

OUU -1250 


Systematic: 

Common- Mode Atmospheric Removal 
Differential-Mode Atmospheric Removal 
Position Offset 
Primary Anisotropics 
Sub-millimeter Galaxies 


+0.06 
-0.06 
+0.01 
-0.04 
+0.12 
-0.00 
+0.05 
-0.05 
+0.20 
-0.20 


+10 

-10 
+50 
-25 
+0 
-0 
+380 
-380 
+650 
-650 


Total: 


o nfi +0.83+0.25 
°- uu -0.83-0.21 


o n n+ 1950 + 755 

ouu -1250-755 



Note. - a The first number is the statistical uncertainty, the second is the 
systematic uncertainty 



-40- 



Table 12. The SuZIE sample of peculiar velocity measurements 



Peculiar Velocity Galactic Coordinates Distance 13 





Redshift 


(km s 1 ) a 


I (deg.) 


b (deg.) 


{h- 1 Mpc) 


A1689 
A2163 
A2261 
A2390 
A1835 
ZW3146 
MS0451 
0016 


0.18 
0.20 
0.22 
0.23 
0.25 
0.29 
0.55 
0.55 


-140011™ 
+1950l|^ 

99^+1650 

-6501?^ 
+750ljgg 
-4050+??° 7 ° 


313.39 
6.75 
55.61 
73.93 
340.38 
239.39 
201.50 
112.55 


61.10 
30.52 
31.86 
-27.83 
60.59 
47.96 
-27.32 
-45.54 


520 
570 
630 
650 
710 
810 
1440 
1440 



Statistical uncertainties only. An additional systematic uncertainty of 
±750kms _1 is assumed for the analysis in §7 

b For a Q m = 0.3, = 0.7 cosmology 



Table 13. Limits to bulk flows in specific directions. 





Direction (Gal.) 


Velocity (km 


s" 1 ) 




Flow name 


1 (deg.) 


b (deg.) 


Best Fit 95% 


conf. limit 


Ref. 


CMB Dipole 


276 


+33 


+60 ± 690 


1410 


1 


ACIF 


343 


+52 


+130 + 510 


1130 


2 


LP10K 


272 


+10 


+40 ± 930 


1860 


3 


SMAC 


260 


+1 


+ 980 


1930 


4 



References. - (1) Kogut et al. (1993), (2) Lauer & Postman (1994), (3) 
Willick (1999), (4) Hudson et al. (1999) 



-41 - 




100 200 300 400 500 

Frequency/ GHz 

Fig. 1. — The frequency dependence of the SZ effect for a cluster with optical depth r = 0.01, 
gas temperature lOkeV and a peculiar velocity of —500km s~ x (towards the observer). The 
thermal SZ spectrum is indicated by the dashed line; the kinematic effect by the dot-dashed 
line. The shaded regions indicate the bands in which SuZIE observes. 
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Fig. 2. — The upper panel shows the layout of pixels in the SuZIE II focal plane. The lower 
panel shows the pattern traced out on the sky by a single column of the array. The cross 
indicates the nominal pointing center of the observation. 
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Fig. 3. — Results of the 2-d likelihood fit of the combined measurements of MS0451 in 
November 1996, 1997, and 2000. The 68.3% and 95.4% confidence regions are shown for 
peak Comptonization and peculiar velocity. 




Fig. 4. — The measured spectrum for each cluster in the SuZIE II sample. In each plot the 
solid line is the best-fit SZ model, the dashed line is the thermal component of the SZ effect 
and the dotted line is the kinematic component of the SZ effect 
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100 200 300 400 100 200 300 400 
v (GHz) v (GHz) 

Fig. 6. — The measured spectra of MS0451 from each of the three observing runs as well as 
the combined spectra using the weighted mean of each spectral point. In each plot the solid 
line is the best-fit SZ model, the dashed line is the thermal component of the SZ effect and 
the dotted line is the kinematic component of the SZ effect. 
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0.7 



Fig. 7. — Left: the stars denote the location, in galactic coordinates, of the clusters observed 
by SuZIE. The direction of the CMB dipole is also shown. Right: the measurements of the 
clusters plotted against redshift. The cross-hatched region shows the range that has been 
probed using optical measurements of peculiar velocities. 
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Fig. 8. — All-sky map showing the upper limit to dipole flow (95% confidence) as a function 
of location on the sky. The locations of the flows listed in Table 13 are also shown. 



